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Abstract 

We apply an improved Taylor expansion method, which is a variational scheme to the Ising model 
in two dimensions. This method enables us to evaluate the free energy and magnetization in strong 
coupling regions from the weak coupling expansion, even in the case of a phase transition. We 
determine the approximate value of the transition point using this scheme. In the presence of an 
external magnetic field, we find both stable and metastable physical states. 



Introduction. More often than not, we encounter situations in which it is difficult to evaluate physical 
quantities by means of standard perturbation methods because, for example, the theory under consid- 
eration does not have any small parameters in which we can expand. As non-perturbative methods, 
variational schemes have been applied to such circumstances with great success. In such methods one 
or more auxiliary parameters are introduced into the model. Optimized perturbation theory [1,2] is a 
systematic improvement of variational methods, formulated on the basis of the "principle of minimal sen- 
sitivity" . [2] In zero and one dimensions, it has been proved that the optimized series converges. [3] This 
method has been applied to, among others, matrix models of superstring theory and their simplified toy 
models. [4-11] It was first reported in Ref. [10] that minimal sensitivity is realized in auxiliary parameter 
space as a plateau (i.e. a region in which physical quantities are stable), and that the appearance of a 
plateau can be regarded as a signal reflecting whether or not the method works. It has also been argued 
that the optimized perturbation theory can be explicitly formulated as an improved Taylor series, which 
is obtained in standard perturbation theory. 

In this article, we attempt to improve our understanding of this method through application to the 
Ising model in two dimensions.^ It is generally believed that physics in an ordered phase cannot be 
analyzed using a naive perturbation theory formulated in a disordered phase, because, in general, a 
perturbation theory yields an ill-defined expansion series when it is applied to a different phase than that 
in which it was constructed. We utilize a mathematical technique to reorganize the series resulting from 
a perturbation theory in a non-trivial manner which is known to yield a good approximation even in the 
case that the expansion is carried out about a point outside the radius of convergence. [10] We show that 
our scheme actually enables us to derive information concerning the ordered phase from a perturbation 
theory formulated in the disordered phase, and vice versa. 

We do not attempt here to calculate the series up to as high an order as possible in order to obtain 
more accurate estimates in the approximation scheme we consider. Rather, our goal is to show that 
physical information can be extracted from relatively low orders of perturbation with this method. In 
actual problems for which the method is intended, such as QED and QCD, no more than first three or 
four terms can be computed in perturbation theory, and from these, we are required to extract physical 
information.^ 

In the following, we first briefly describe how the optimized perturbation method can be understood 

as an improved Taylor expansion (ITE). Then we can apply it to the Ising model in two dimensions 
to determine whether the method works in a system that exhibits a phase transition. By applying the 
ITE to this model, we show that the ITE provides a phase which is invisible from another phase in 
the perturbation expansion. We also investigate within this model the extent to which we are able to 
understand the critical behavior from the ITE analysis. 

Improved Taylor expansion. Suppose we have some quantity F{X, m) which is an implicit function 
of the parameters m and A, and we would like to evaluate it at some definite values of these parameters. In 
order to treat F{X, m) using perturbation techniques, we regard one of the parameters, A, as an expansion 
parameter. However, at this point we do not assume that A is necessarily small. Also, in general, the 
parameter m might respect a set of parameters, but we will not consider such cases in this article. Then, 
one might perform a formal perturbative expansion with respect to A up to some finite order, say, A''th 
order, which depends on the nature of the problem at hand (or our ability to carry out the necessary 
calculations). From the A'' -|- 1 coefficients appearing in the expansion, we obtain a sequence of finite 
series Fq,Fi, - ■ ■ , Fjv in the form 



^In Ref. [13] there appears a treatment of the application of the improved Taylor expansion method to the Ising model. 
Our arguments here extend and clarify the work given there. 

^One attempt to apply this method to non-abelian Yang-Mills theory can be found in Ref. [14]. 
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It is not necessarily the case that each such series is convergent, because the expansion parameter A might 
be large (outside the radius of convergence), or even for small A, it could turn out to diverge at some 
order if it is an asymptotic series. Our goal here is to improve the convergence of this series. For this 
purpose, we first introduce a formal expansion parameter g and an auxiliary parameter mo, and then 
make the substitutions 

A gA, 

m mo + g{m — toq) , (2) 
in the arguments of each series Fi(X,m). Then, it is easily seen that 

i i 

J2 ^"/"(™) = E 9">^".fn{mo + g{m ~ mo)) 



n=0 n=0 



9=1 



n=0 \k=0 



(3) 

9=1 



where the notation \g_i means that we substituted 5 = 1 at the end. It might not be obvious from the 
above expression that the last expression of Q is independent of the parameter mo as a whole, but in 
fact they are. However, if we drop the 0{g^~^^) terms, the coefficient functions become dependent on mo. 
Following these steps, we obtain a new sequence of series which we call the "improved functions" : 

Fr''(A.m;m„).i:A"(g ''"';"°''"'t',r' ). (^^o,.^.iv) (4) 

n=0 \k=0 ' / 

According to the principle of minimal sensitivity, [2] the improved functions best approximate the exact 
value when they are evaluated in regions where they are least sensitive to the auxiliary parameter mo. 
We refer to the regions in which the minimal sensitivity condition is realized as plateaux. [10] 

It is a rather difficult problem to identify a plateau, especially when there are many auxiliary param- 
eters.^ Indeed, we have no quantitative nor mathematically rigorous definition of a plateau. However, in 
the case of a single auxiliary parameter, one may plot graphs of the improved series as functions of the 
parameters. Then, examining these graphs, we are able, at least qualitatively, to determine whether or 
not a plateau exists. 



Ising model. The Ising model [16] is one of the simplest models that exhibit phase transitions. Here, 
we test the ITE using the Ising model on a two-dimensional square lattice with a uniform magnetic field 
h. The total number of sites is denoted by V, which defines the volume of the system. The Hamiltonian 
is given by 

= - J ^ O-jCTj - /l ^ (Ti , (5) 

where ct^ located at the site labeled by the index i takes one of the values ±1, and there are interactions 
only between nearest neighbors, with the strength J (> 0). Pairs of neighboring sites are denoted as (ij). 
The partition function is defined as 

Z(J,;i)^^e-^, (6) 
W} 

where we normalize the coupling and external field strength in units of fcsT. The sum is taken over all 
possible configurations. 

''In Refs. [10-12] the accumulation of extremal points in the parameter space is used as a criterion for the identification 
of a plateau. An efficient method to detect a plateau in a multi-dimensional parameter space is presented in Ref. [6], where 
the histogram method is proposed. Further discussion will bo given in a forthcoming paper. [15] 
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The exact free energy of the two-dimensional Ising model was obtained by Onsager: [17] 



/exact (^) =ln 
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~8n^J J l^[(l + ^^)^-2A;(l-fc2)(cosa; + cosy)] , (7) 
where k = tanh(J). Also, the exact expression of the magnetization [18] is 

m(J) = [l-sinh-4 (2 J)] (8) 

In the following, we first apply the ITE to the weak coupling expansion and then to the strong external 
magnetic field expansion. 

Weak coupling expansion. In this section, we employ the weak coupling expansion* and construct 
the following series: 

oo 

Z(J,/i) = 5^J"^„(/i). (9) 

n=0 

Calculating 2„(/i) up to fourth order in J we have the following: 

2o(/i) = (2cosh(/i))^, (10a) 
zi{h) = zo2Vtanh^{h), (10b) 
Z2{h) = ^2V [{2V - 7) tanh*(/i) + 6tanh^(/i) + l] , (10c) 

Z3{h) = ^2V [(4^2 _ 42V + 116) tanh^(/i) + (36F - 168) tanh*(/i) 

+{6V + 52) tanh^(/i)] , (lOd) 
Z4{h) = ^2V [{8V^ - 168^2 + i222y - 3126) tanh^(/i) + (144^^ _ i848V 

+6312) tanh^(/i) + (241/^ + 548^ - 3748) tanh''(/i) + (72^ + 552) ta.nh^{h) 

+{6V + 10)] . (lOe) 

Note here we are ignoring the effect of the boundaries, which is irrelevant in the large volume (thermo- 
dynamic) limit. 

The free energy density at jth order, fj{J, h), is obtained as a series expansion up to order j of the 
logarithm of the partition function: 



(11) 



Here, the notation 



indicates that all terms of order j + 1 or higher in J are omitted. Because we have 

calculated the partition function up to fourth order in J, we can easily obtain the fourth-order free energy 
density as a function of J and h: 

U{J,h) = -ln(2cosh(/i)) -2tanh^(/i)J- {1 -F6tanh2(/i) - 7tanh''(/i)} 

y tanh2(/i) - 56tanh''(/i) + ^ tanh®(/i)| ^^2) 

^ QS7 521 1 

- + 46tanh^(/i) - — tanh*(/i) + 526tanh^(/i) - — tanh^(/i) \ J* . 

Do ^ J 



This series can be understood as the high temperature expansion in the case of zero external magnetic field, h = 0. 
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This function is independent of V, as should be the case. 

The magnetization density is obtained by differentiating the free energy density with respect to the 
external magnetic field h: 

mj{J,h) = — . (13) 

This leads to 



mi{J,h) = tanh(/i) 



1 + 4{1 - tanh(/i)^} J + 4{3 - lOtanhihf + 7tan]i{hf}J^ 



tanh(/i)2 + 456 tanh(/i)'* - 232 ta.nh{h)^ | (14) 



104 776 
"3 3~ 

92 _ tanh(/i)2 + tanh(/i)4 - 5240 tanh(/i)*^ + 2084 tanh(/i)** I J"* 

In the following, we refer to the free energy density and magnetization density as simply the free energy 
and magnetization. 



Applying the ITE. Having thus obtained the basic ingredients through perturbative calculations 
about J = 0, we now apply the ITE to the resulting series in order to obtain improved information, 
especially at large J (low temperature). For this purpose, we introduce a formal expansion parameter g 
and an auxiliary parameter /iq, and then make the substitutions 

J ^ gJ , 

h^ho + gih-ho) (15) 

in the expressions for the free energy, (|12|l . and magnetization, H14|) . After carrying out a Taylor expansion 
about g = up to the specified order and then setting g = 1, we obtain improved functions. In the present 
case, we have calculated the free energy and the magnetization up to fourth order in g. We have sequences 
containing five improved functions of the free energy, f™"^ {J , h] ho) , and magnetization, m™'^{J,h;hQ), 
which contain the fictitious parameter Hq. We list the explicit forms of these improved functions in 
Appendix 1X1 

In order to investigate the spontaneous magnetization, we set the external magnetic field, h, to zero 
in the improved functions. The results obtained using the ITE can then be compared with the exact 
values, {Tj) and (|SJ|, to test the validity of the method. 

Before discussing the improved functions in detail, let us first comment on the behavior of the functions 
in the case — and — s- ±oo. For Hq — 0, the improved free energy at ith order coincides with the 
original one at the same order as a function of J; e.g., we have 

/r''( J, 0; 0) = - In 2 -J2-^J4 

= /4(J,0). (16) 

For the asymptotic behavior in the limits ho — > ±oo, all of the improved free energies become close to 
—2 J (except /g™^ which diverges): 

fl"'P{J,h-ho^±oo)^±h-2J, (i = 1,2,3 and 4) (17) 

where we have restored h for later convenience. This shows that the improved free energies are insensitive 
to the value of the parameter ho for large \ho\, and thus, two flat regions are formed. However, we do 
not regard those asymptotic flat regions as plateaux, because their forms are independent of the order 
of perturbation; we consider a flat region to be a plateau only if its form changes as the order increases. 
This also implies that in order to determine whether a given region is a plateau, we have to consider the 
whole sequence of improved series, rather than only the functions of highest order. 
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Figure 1: (a) Dependences of the fourth-order improved free energies, (ho), with respect to the 



auxihary parameter ho for J ~ 0.01, 0.1, 0.2, 0.3, 0.4, and 0.5. Jh) Dependences of the improved free 
energies fl"^^{ho) for J = 0.6 from zeroth order to fourth order. The horizontal hues represent the exact 
values of the free energies at the respective values of J. 



With the above points in mind, let us study the behavior of the improved functions in detail. Figure 
1(a) exhibits the fourth-order improved free energy J™^ as a function of ho for J = 0.01, 0.1, 0.2, 0.3, 0.4, 



and 0.5. As we vary J, there is a clear evolution of the shape of the improved free energy. This reflects the 
fact that there are two phases, corresponding to the small J (disordered) phase and the large J (ordered) 
phase in the Ising model. At small values of J (i.e. in the high temperature region), the set of improved 
free energies form a clear plateau around ho = 0, at which the values of the improved free energies are 
good approximations of the exact value. This is not surprising, because the original free energy was 
obtained through a perturbative expansion in small J, and it already provides a good approximation of 
the exact value. As we have seen above, at ho — 0, which is on the plateau, the fourth-order improved 
free energy has the same form as the original. 

Now, let us consider the case of relatively large J (i.e. the low temperature region). In this case, it is 
seen that plateaux develop in the regions near ho ^ ±2 as the order of the improved free energy increases. 
Thus, we can evaluate the approximate value of the free energy on these plateaux, which is found to be 
close to the exact value expressed by the horizontal line [Fig. l(b)| . We would like to emphasize that we 
have obtained the relatively large J (ordered phase) physics from the data of small J (disordered phase) 
expansions. 

Although the emergence of plateaux at non-zero ho in the improved free energies may signal the 
occurrence of spontaneous magnetization, we have to investigate the magnetization itself to determine 
the amount by which the system magnetizes spontaneously. In Figs. |21 and |31 respectively, we plot the 
improved magnetizations from zeroth to fourth order for J = 0.1 and J — 0.6. In Fig. [3 it is seen that 
the slopes become gentle and a plateau emerges near the origin as the order increases. Contrastingly, 



Fig. 3(a) shows that the slopes at the origin become steeper as the order increases, whereas two plateaux 
develop near ho — ±2. This figure also provides an approximate value of the magnetization that is close 



to the exact value, represented by the horizontal line in Fig. 3(b) 



The two phases are characterized by spontaneous magnetization. If the plateau corresponds to the 
ground state realized at a particular value of J, the critical point can be recognized by the disappearance 
of the plateau formed at and near the origin. As mentioned above, the slope at the origin becomes 
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Figure 2: Dependences of the improved magnetizations m""'^{h()) with respect to the auxihary parameter 
ho at J — 0.1. 
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Figure 3: (a) Dependences of the improved magnetizations m™^(/io) with respect to the auxihary pa- 
rameter ho for J = 0.6. (b) An enlarged plot near the exact value, which is represented by the horizontal 
line. 



steeper or gentler as we increase the order, depending on the value of J. Thus, at each order of approx- 
imation, the critical point is determined by the condition that the difference between the derivatives of 
the magnetizations at nth and (n — l)th orders be zero at the origin. Thus, we find^ 
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( J(i) = 0.25) 
( Jp) 0.33333) 
(4^) 0.34615) 
(4^^ ^ 0.37681) 



(18a) 
(18b) 
(18c) 
(18d) 



It is seen that as we increase the order, the approximate value of the critical point moves closer to the 
exact value, Jc = 5 ln(l + \/2) = 0.440 • • • . 

It is interesting that the quantity on the left-hand side of (|18|l is equal to the second derivative of the 
improved free energy with respect to Hq at the origin: 



2 fimp 



d'f, 



dhl 



ho=0 



dm. 



imp 



dho 



o imp 

dho 



(19) 



J ho=0 



Thus, the phase transition point can be understood in terms of the improved free energy in the following 
way. The plateau at the origin in the disordered phase is a local minimum, and as we increase J, it 
becomes a local maximum, and new minima appear. The critical point is the phase transition point. 

It is also interesting to see what happens when we turn on an external magnetic field. For this 
purpose, we substitute h = 0.05, as an example, in the improved functions. Turning on the external 
field causes the heights of two degenerate plateaux to become different [Fig. 4(a) . The plateau whose 



height increases corresponds to a metastable physical state (vacuum) , and the other plateau corresponds 
to a stable physical state. Hence, the ITE reveals not only a stable vacuum but also a metastable one. 
However, as we increase the strength of the external field, the metastable plateau disappears, while the 
stable plateau retains its flat shape. 



®The critical temperature at first order (Jc^^ = 0.25) is the same as that calculated in the mean field approximation. 



7 













ordorl 




1 '■ 






order 3 




order 4 




/ 
/ 


\ exact(h=0) 



























order 2 
" order 3 " 

<a°[| hLol 



ho 



(a) 



(b) 



Figure 4: (a) Dependences of the improved free energies f^ (ha) with respect to the auxihary parameter 



ho in the case of a weak external field {h = 0.05) and J = 0.6. Jh) An enlarged plot. The horizontal line 
represents the exact value in the case h = 0. 

Strong field expansion. Let us now proceed to the strong field expansion^ and construct a series in 
terms oi H = exp(— 2/i): 



Z(J,/i) = ^i?"z„(/i). 



We calculate z„ up to fourth order in H: 



^V{2J+h) 



zo{h) 



Z2{h) = ZqV 

z^ih) = z^V 
Zi{h) = zoV 



A* + 18A^ + (8y - 85)A'^ + {V^ - 21V + 96)A^ 



+ - {V^ - 30V"2 ^ 323^ + 714) 
where A = exp(— 4J). The free energy density of jth order is defined as 



■lln {^H-z^{h)^ 



(20) 

(21a) 
(21b) 

(21c) 
(21d) 

(21e) 
(22) 



^The strong field expansion is closely related to the strong coupling expansion. Actually, in the strong coupling expansion, 
a small external field is introduced to split the degeneracy and to decouple one of the sectors in the thermodynamic limit. 
The remaining sector is essentially that obtained from the strong field expansion. [19] 
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Figure 5: (a) Dependences of the fourth-order improved free energies fl"^^{H = 1, J; Ao). (b) Those of the 



improved magnetization rn^^{H = 1, J; Ag) in the strong field expansion for J = 0.3, 0.4, 0.5, 0.6, and 0.7. 



We obtain the free energy of fourth order as 

/4(J, h) = -(2J + h)- X^H - (^2A3 - ^a4^ _ (^QX^ - 16A5 + ^A^ ) 

121 \ 

A* + 18A^ - 85A^ + 96A^ ~ — A* | 



(23) 



and the magnetization as 



m4( J, h) = l- 2X^H + 4 (^-2A3 + ^A*^ + 6 (^-6X^ + 16X^ - y A^ ) 

.A4 - iSA^ + 85A6 - 96A^ + — A« ] . 
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(24) 



In order to obtain improved functions, we make the substitutions 



A -> Ao + g{X - Ao) 



(25) 



in the above expressions for the free energy and magnetization. Expanding about g up to fourth 
order and setting 5 = 1 at the end, we obtain the improved free energy /™^(i/, J; Ao) and magnetization 
m™^^{H, J; Ao). We set H = 1 because we are interested in the case in which the external magnetic field is 
turned off. Here we present the fourth-order improved free energy /™^(ff = 1, J; Ao) and magnetization 



\H = 1, J; Ao) for various J in Figs. 5(a 



and 



5(b) 



It is curious that there is no clear transition in 



the shape of the function as we vary the strength of J, in contrast to the previous case. Nonetheless, an 
indistinct signal can be seen that the plateau located near Aq = moves to Aq > below a certain value 
of J, which may imply the occurrence of a phase transition. 



Conclusion. We applied the improved Taylor expansion (ITE) to the Ising model in two dimensions 
and found that it is effective for extracting quantities in strong coupling regions from the weak coupling 
expansion, even though these two regions are separated by a phase transition. We obtained a good 
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approximation for the observable in the ordered phase from a perturbation theory formulated in the 
disordered phase. We also determined approximate values of the critical point at each order of the 
approximation by closely examining improved functions of the magnetization, in particular at the origin 
of the auxiliary parameter space. Furthermore, turning on a small external field the ITE reveals the 
existence of not only a stable vacuum but also a metastable one. We also applied the ITE to the strong 
field expansion and found that it also works well for evaluating the free energy and magnetization in weak 
coupling regions, oven though there is only a weak signal indicating the phase transition, which can be 
seen clearly in the weak coupling expansion. 

We note that the ITE at first order is identical to the ordinary mean field approximation when we adopt 
a criterion for the identification of the plateau using the stationary points of the improved free energy. The 
condition that the derivative of the improved free energy with respect to the auxiliary parameter be zero 
should here play the role of the self-consistency condition in the mean field approximation. Extensions of 
the mean field approximation have been explored, e.g., in the Bethe approximation, in such a way that 
the quantum degrees of freedom of neighboring sites are taken into account step by step. This approach 
incorporates the higher order contributions successively, while the physical picture is apparent. Our ITE 
scheme also provides a systematic extension in a simple manner based on the variational method in a 
purely mathematical sense. It turns out, however, that the ITE is different from Bethe approximation 
already at second order, and the physical meaning of the extension is not obvious. It would be interesting 
to clarify the role of the plateau and its relation to the self-consistency condition in the mean field 
approximation. It would also be interesting to evaluate various critical exponents within the ITE scheme 
as a future work. 
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stage of this work. T. M. also thanks T. Koretsune for valuable discussions. T. M. and Y. S. are supported 
by the Special Postdoctoral Researchers Program at RIKEN. 
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A Improved Functions 

Here we list the improved free energy and magnetization up to fourth order in the weak coupling expan- 
sion. In the following, we employ the quantity ko = tanh(/io)- 

ft" {J, h; ho) = - ln(2 cosh(/io)), (26a) 
fr^iJ, h; ho) = -{ln(2 cosh(/io)) + ko{h - ho)} - 2fc^ J, (26b) 

fr^iJ, h; ho) = - |ln(2cosh(/io)) + ko{h- ho) + ^(1 - kl){h - /io)'| 

- {(4fco - 4fc3)(ft - ho) + 2fc2} J _ {1 + 6fc2 - 7fc*} J2, (26c) 

/r^(J, h; ho) = - |ln(2 cosh(/io)) + ko{h - ho) + ^(1 - kl){h - ho)^ 

-^(fco - kl){h - ho)A - {2^0 + (4fco - 4fc^)(/i - ho) 



3 

+ (2 - 8kl + 6k^)(h - ho)^}J - {1 + 6fc2 - 7k^ 

f ^2 nfi 

+ (12A;o - AOkl + 28k^){h - ho)}J^ - <^ — fcg - 56fc^ + — fcg } J^ (26d) 



/r^(J, h; ho) = - |ln(2cosh(/io)) + fco(/i - ho) + ^(1 - fc^)(/i - ho)'' 
— {ko - kl){h - ho)' + (-^ + \kl - ^fco') {h - hoY 

2fc2 + (4fco - Akl){h - ho) + (2 - 8fc^ + 6fc^)(/i - /lo)^ 
+ (-yfco + yfco' - 8fc^) {h - /io)'| J - {1 + 6fc2 - 



+ (12fco - 40fc^ + 2S,kl){h - ho) + (6 - 66fc^ + 130fc^ - 70k^){h - /lo)^} 

-{I 

6 ■ 3 ■ ^ 2 



^ ^^fc2 - 56fc4 + li^fcS ^ 1^104^^ _ 776^3 ^ ^^^^^ _ ^g^^r^ _ ^^^j 
5 + 46fc2 - + 526A;6 - ^kl] J^ (26e) 
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rn^^{J,h; ho) = ko, 

mr^(J, h; ho) = {ko + (1 - - ho)} + {Ako - Ak^jJ, 

mf ^(J, h; ho) = {ko + (1 - kl){h - ho) + {-ko + kl){h - ho)^} 

+ {4/00 - 4fc^ + (4 - 16kl + 12k^){h - ho)} J + {12ko - 40fc^ + 2Sk^}J^, 

mt"P{J, h; ho) = |/eo + (1 - ki){h - ho) + {-ko + kl){h - ho)^ 

+ (-1 + - k^) {h - ho)'} + {iko - 4fc3 + (4 - 16fc2 + l2k^){h - h 



.33 

+ {-l6ko + 40fci^ - 2Ak^){h - ho)^}J 

+ {12fco - 40fc3 + 28fc^ + (12 - 132fc^ + 260fc^ - UOk^){h - ho)}J'^ 
+ I ^fco - ^-f-kl + 456fco^ - 232/e^l j\ 

r^r{J, h; ho) = |a-o + (1 - kl){h - ho) + {-ko + kl){h - ho)^ 

+ (-1 + - kt^ {h - ho)' + - Ik'o + kl^ {h - ho)'^ 

+ I^Ako - 4:kl + (4 - 16kl + 12k^){h - ho) + {-Uko + 40k^ - 24:k^){h ■ 

y + ^k'o - 80kt + AOk{^ {h - ho)'^ J 

+ {12fco - 40fc^ + 28fc^ + (12 - 132fc^ + 260k^ - UOk^){h - ho) 
+ (-132fco + 652fc^ - MOk^ + 420fcJ)(/i - ho)^}J'^ 

^ko - + 456fc5 - 232^0^ 



104 2432 



3 3 



fcg + 3056A;^ - 3904fcg + 1624fc^^ {h-ho)^J' 



+ j^92ko - ^k'o + ^^k^o - 5240fco" + 2084fc9| j\ 
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